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ABSTRACT 



We analytically compute the long-term variations of all the standard six Ke- 
plerian orbital elements a, e, /, Q, w,Ai of a test particle orbiting a central body 
acted upon by an incident monochromatic plane gravitational wave. We assume 
that the characteristic size of the perturbed two-body system is much smaller 
than the wavelength of the wave. Moreover, we also suppose that the wave's fre- 
quency Ug is much smaller than the particle's orbital one n. We make neither a 
priori assumptions about the direction of the wavevector k nor on the orbital ge- 
ometry of the planet. We find that, while the semi-major axis a is left unaffected, 
the eccentricity e, the inclination /, the longitude of the ascending node Q, the 
longitude of pericenter w and the mean anomaly Ai undergo non- vanishing long- 
term changes of the form d^/dt = (S(/Cjj; e, /, f2, w), \& = e, I,Q,zu, Ai where 
JCij, i,j = 1,2,3 are the coefficients of the tidal matrix K. They are not secu- 
lar trends because of the slow modulation introduced by K and by the orbital 
elements themselves. Our calculation holds, in general, for any gravitationally 
bound two-body system whose characteristic frequency n is much larger than the 
frequency Ug of the external wave. It is also valid for a generic perturbation of 
tidal type with constant coefficients over timescales of the order of the orbital 
period of the perturbed particle. 

Subject headings: gravitation; gravitational waves; celestial mechanics 
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1. Introduction 



Gravitational waves are a key theoretical prediction of the Einstein's theory of 
gravitation. Indeed, since the general theory of relativity relies upon the Lorentz invariance, 
which carries with it the concept of a limiting sp eed for interactions, the existence of 
gravitational waves is a natural consequence of it ( Einstein 1916 , 19181 ). Apart from the 
indirect evidence of their existence infe rred from the orbital decay rate of the binary pulsar 
PSR B19 13+16 (iHulse fc Tavlorl Il975l). a di rect measurement of them is still lacking. 
See, e.g.. ICerdoniol fl2003[ ): iGiazottd (120081 ): iFairhurst et al.l (120111 ) for recent reviews of 
the status of gravitational wave detection. The consequences of detecting gravitational 
waves for physics, astrophysics , and cosmology would be certainly remarkable; see, e.g., 
Sathyaprakash fc Schutd ( 20091 ). The role of a direct detection of the gravitational waves 



for general rel ativity and extended theories of gravity has been recently pointed out by 
Cordal ( l2009al Jbl). Let us recall that the gravitational wave spectrum spans a wavelength 
range of over 18 orders of magn itude, encompassing a very broad range of physics and 
astrophysics (jPrince et al.ll2010l ). The frequen cies in the range 10^ — 1 0^ Hz are the targe ts 
of several ground-based detectors like, LIGO ( Abramovici et al. 1992 : Abbott et al. 2009h. 



VIRGO (Bradaschia et al. 1990: Vinet 



Takahashi et al. 2004: Aral et al. 2008 



)■ GEO (iDanzmann et al.lll992l : iGrote et al.ll2008f ) 



19931: lAcernese et al.lboo J. TA MA (lTsubono[ 



1995 



etc. The space-based LISA mission (jShaddockll2008l ) aims to detect gravitational waves in 
the frequency range 10"^ — 1 Hz, whi le accurate tim i ng measurements of pulsars rn ay detect 
signals in the range 10^^ - lO"'^ Hz (lKopeikirJll997t lHandll2010l : Ijenet et al.lboilh . 



In this paper, we will analytically work out the long-term variations of all the six 
standard Keplerian orbital elements of a solar system planet due to the action of an 
externally incident monochromatic plane gravitational wave in the green-recQ part of the 
spectrum, i.e. with frequency Ug <C 10~^ — 10~^° Hz. Such kinds of gravitational waves are 
important since they carry information about how galaxies and blac k holes co-evolved over 
the history of the Universe (jPlowman et al.ll2010l : ISesana et al.ll201ll ). the early universe and 
related exotic physical processes like, e.g., inflation and cosrnic strings, and possible phy sics 



beyond the standard model of particles and fields (IVilenkinI Il98ll : iRubakov et al.l Il982 



Fabbri fc Pollocklll983t lHoganlll986l : lGrishchuklll993l ) . Concerning the possible existe nce of a 



(1961 


): 


Wheeler 


Mashhoon et al. 



(1962) 
(Il98lf ) 



Zel'dovich fc Novikovl ( 1l975l ) ; ICarrI (jl980[ ) , and the discussion in 



The effects of incident gravitational waves on the orbital motion of 
gravitationally bound systems were considered by several authors with a variety of different 



themselves 


(Bertotti 


19731: 


Rudenkol 


1975 


Mashhoon 


1978 


Mashhoon 


L979|: 


Turnerl 


19791: 


Grishchuk & Polnarev 


1980|: 



See Figure 1 of Prince et al. (l2010f ) 
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1982a 


b 


1996a 


b 



Nelson & Chau 


1982; 


Ismaiel & Saad 


2011) 



Ivashchenkd 119871 : iKochkin fc Sbytovi Il987l : IChicone et al. 



stochastic background of gravitationa l waves of wavelengths much larger than about 1 au 



was first suggested by iBertottil (119731 ). We wish to remark that the calculation presented 
here are, actually, of quite general interest. On the one hand, as far as the search for 
gravitational waves is concerned, our results may b e extended also to spac ecraft-based 



missions in the solar system like ASTROD-GW (INil 120091 : Men et all 120101). const ituted 



of a number of probes in wide heliocentric orbits, and LISA ( iPovoleri fc Kembld 12006 



Xia et al.ll2010l ) envisaging the use of three spacecrafts orbiting the Sun at 1 au distance in a 
quasi-equilateral triangle formation 20 deg behind the Earth. It is also the case of recalling 
that the possibilities offered by Doppler tracking of interplanetary drag-free spacecraft to 



detect cosmological gravitat ional waves wit 



au w ere studied in the past flAndersod 11971 



1 wavelength of the order of, or larger than 1 



Mashhoon fc Grishchuk 1980: Bertotti fc Carr 



1980l ): for a recent review, see lArmstrongl (120061 ) and references therein. Moreover, our 



calculations are va lid, in principle, also for other natural systems like extrasolar planet^ 
(ITorres et al.ll2008l ) many of which have orbital frequencies of the order of 10~^ Hz, orbiting 
their parent stars at distances as small as ^ 0.01 au. In such cases, our findings are 
technically valid for waves with higher frequencies with respect the solar system ones: 
indeed, they might be as hi gh as abouijl 1 0~^ co rresponding to the green-light blue part of 
the spectrum in Figure 1 of iPrince et al.l (120101 ). On the other hand, our results are not 



necessarily limited to the very-low frequency waves case, being valid for any tidal force with 
constant (over particle's characteristic timescales) matrix coefficients as well, independently 
of its physical origin. 

The plan of the paper is as follows. In Section [2] we will discuss the analytical form 
of the acceleration experienced by the test particle due to an incident monochromatic 
plane gravitational wave traveling along a generic direction. We will also consider the 
simplified cases of a gravitational wave moving along the z axis, as it is a choice widely 
adopted in literature, and in the reference {x, y\ plane. In Section [3] we will analytically 
work out the long-period changes occurring in the particle's orbital motion in the case of 
wave's frequencies much smaller than the orbital ones. We will not make any simplifying 
assumptions about the inclination and the eccentricity of the orbit. In Section H] we briefiy 
review some of the approaches followed in literature. Section [5] is devoted to the conclusions. 



^See |http://exoplanet.eu/ on the WEB. 

■^However, for waves with such relatively high frequencies other facilities like LISA would 
be available: if they will finally become operative at the expected level of accuracy, they 
would likely surpass the possibihties offered by the extrasolar planets. 
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The acceleration imparted on an orbiting planet by a passing 
monochromatic plane gravitational wave 



The action of an incoming gravitational wave on a planet of the solar system is of tidal 
type with respect to a suitabl y constructed local inertial frame, represen ted by a Fermi 
coordinate system {x,y,z,t} ( Fermilll922l : iLevi-Civital Il926l : ISyngelll960[ ). whose origin is 
located at the Sun's position. In general, a tidal acceleration A experienced by a slowly 
moving test particle due to an external curved spacetime metric can be written in terms of 
the "electric" components R\jo, i,j = 1,2,3 of the local Riemann curvature tensor R as 
(jPiranilEoseh 



— 



1,2,3, 



where we introduced the coefficients K, 



Levi-Civita 


1926; 


Svnee 


1926) 



3 of the tidal matrix K. In this case, the 



In general, other non-linear terms involving also the components of the test particle's 
velocity along wi th other coinpone r its of the Riema i in tensor occur from the generaliz ed 
Jacobi equation f lHo dgkinsonl 1 1 9 721 : lMashhoonlll975l . 119771 : IChicone fc Mashhoonl[2002l ): we 



will neglect them because of their smallness. 



In the linearized weak- field and slow- motion approximation (jPadmanabhanI 120101 ). and 
for the case of a propagating gravitational wave, one ha^ 



R 



1 d'^h, 



2 df^ 



,^,j = 1,2,3. 



(2) 



Recall that, in general, the symmetric tensor h (jPadmanabhanI |20 1 Ol ) 



0,1,2,3, 



(3) 



where h^^ , /i, z/ = 0, 1, 2, 3 are small {\h^^\ <^ 1, /i, = 0, 1, 2, 3) corrections of the metric 
tensor 77'^'^, /i, z/ = 0, 1, 2, 3 of the fiat Minkowskian spacetime, and h is the trace of h, is a 
solution of the linear wave equation in the vacuum 



0, fi,u = 0,l,2,3, 



(4) 



with the gauge condition 



■flU 



dh 
dxf^ 



0, z/ = 0,1,2,3. 



(5) 



^In general, also the second derivatives of /iqo and h 



Oi, 



1, 2, 3 are present in R\,n, i, j 



1,2,3, but in the case of a gravitational wave they do not propagate (jPadmanabhanI |2010[ ) 
because of gauge freedom (transverse traceless gauge, or TT gauge). 
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In the TT gauge, it turns out that h = h. In the case of a propagating plane wave, of the 
form@ 

h'"' = Re {x'"" exp [i {ugt - fc • r)]} , /z, i/ = 0, 1, 2, 3, (6) 

the fact that both /iqo and hoi, i = 1,2,3 vanish in the TT gauge imphes that the 
tidal matrix K is not only symmetric, but also traceless. Thus, it has five independent 
components, so that the acceleration A of eq. ([1]) felt by a test particle becomes, quite 
generally. 



(7) 



In fact, the condition that the force exerted on the particle by the wave is orthogonal to its 
direction of propagation given by k yields three further constraints 



'^xxf^x ' i^xy'^y ' '^xz'^z 



^xy^x ~1~ ^yy^y ^yz^z 0; 



^ ^xz^x ~l~ ^yz^y {]^xx ~l~ ^yy^ ^z 0; 



(8) 



so that the independent components of ICij, i,j = 1, 2, 3 are just two. They have dimensions 
of T~^ and, in the case of eq. ([6]), are harmonic functions proportional to u^. Notice that 
eq. ([8]) comes just from eq. ([5]): indeed, for eq. ([6]) and in the TT gauge eq. ([5]) reduces to 



(9) 



which corresponds to eq. ([8]). 

For a wave not traveling in the reference {x, y} plane, i.e. for kz 7^ 0, eq. (IE]) yields 



yy 



1C..X ( - 2IC„, ( Mh) 



yz 



IC^y < f I- 



(10) 



2kl 



\ky+k^j 



Hn eq. (jH]) x'^'^^ fi,^ = 0,1,2,3 represent the gravitational radiation field. 
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Thus, we pose 



^1 ^XXI 

h2 = A^T 



(11) 



for the two independent polarizations of the gravitational wave. 



When the wave propagates in the reference {x, y} plane, i.e. for = 0, the second 
equation in eq. fllOp becomes 



^yy 



2/C,, 



xy 



(12) 



while from eq. (|H]) it turns out 



/C 



xy 



-IC ^ 



(13) 



/C 



yz 



Thus, in this case, we define 



hi = 



/C.T 



(14) 



Notice that both eq. f[T^ and eq. f[T^ do not hold for a wave traveling along the reference 



X axis, i.e. for ky = 0. 



Finally, let us notice that when k = {±1, 0, 0}, eq. (j8]) tells us that 



so that it can be posed 



^xx ^xy ^xz 0, 



^2 = !^yz- 



(15) 



(16) 



In order to make contact with realistic situations occurring in typical solar system data 
analyses, we remark that the unit vector k can be written, in general, as 



kx = cos /3 cos A, 



< ky = cos 13 sin A, 



(17) 



, kz = sin /3 
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in terms of the ecliptic latitude /3 and longitude A, which are to be considered as unknown. 
Explicit expressions of the tid al matrix coefficients for a generic wave's incidence can be 
found in lChicone et al.l (jl996al ). In addition to the amplitudes of the two independent wave's 
polarizations and of their mutual constant phase difference, they contains f3 and A through 

the angle (IChicone et al.lll996a| ). It is defined from k • N = cosG and k x N = sin0, 

where N is the unit vector directed along the test particle's out-of-plane direction (see eq. 
( 13T|) below) coinciding with the direction of the orbital angular momentum. 

Inserting eq. ( TTUj) . with eq. (fT7|) . in eq. (j7]) yie 



'Ax = hi{x — z cot /3 cos X) + h2{y — z cot /3 sin A) 



A, 



A. 



j sin^ /3+cos^ /3 cos^ A . , 



cot l3 sin A ( sin cos A I 

sin"^ /3+cos^ /3 sin"^ A 



+ 



+ h2 



. cos-^ /3 sin 2A 

~ y I sin^ /3+cos2 13 sin^ A 



z cot /3 COS A f 

\ sin'^ ^+cos^ 13 sm^ A J 



hi < —X cot /3 COS A + y 



cot (3 sin A (sin-^ /3+cos^ (3 cos^ A) 
sin^ /3+cos2 (3 sin^ A 



I I cos^ f3 cos 2A 
~'~ ' sin^ /3+cos2 /3 sin^ A 



+ h2 



-X cot /3 sin A - y cot /3 cos A ( ^7^°^^ ^ ^ ) + ^ ( . /^I'^f/^.. 

^ » \ sin-^ /3+cos^ /3 sin*^ A y \ sin*^ /3+cos"' /3 sin^ A 



In the specific case of a plane wave propagating along the reference z axis, i.e. for /3 = ±90 
deg, eq. ^ and eq. (fTOj) yield 



r ^ 



/c 



(19) 



0. 



Thus, eq. (fT8l) reduc es to the well known result ( jNelson fc Chaul Il982[ llvashchenkd Il987 
Ismaiel fc Saadll201ll ) 

= hix + /i2y, 



A, 



h2X - hiu, 
0. 



(20) 



^Notice that eq. ( !T8|) has a singularity for /3 = 0, i.e. for a wave traveling in the reference 
{x, y} plane. 
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For 13 



0, i.e. for kz = 0, eq. ([7]) reduces tc0 

'Ax = hi{x — y cot A) + h2Z, 

Ay = — cot A [hi {x — y cot A) + h2z] 

A. = hz {'^) + h2{x-y cot X) 



(21) 



because of eq. f lT2|) . eq. f lT3|) and eq. f lTTj) . Notice that eq. f l2T|) is not singular for 
= {0, ±1, 0}, and A^ 0,Ay = 0, 7^ 0, as it is expected for a plane wave. 

Finally, for a wave traveling along the reference x axis eq. (^^, with eq. ( IT^ and eq. 
f lTB]) . becomes 

r A. = 0, 



A,, 



h2y - hiz. 



(22) 



In Section |3] we will analytically work out the effects of eq. ([7]) on the trajectory of 
a test particle orbiting a central body with gravitational parameter GM, where G is the 
Newtonian gravitational constant and M is its mass, located at the origin of the Fermi 
frame traversed by a monochromatic plane gravitational wave. We will also consider the 
particular case of of eq. ( 12U]) (fc = {0, 0, ±1}), widely treated in literature. 



3. The long-term variations of the Keplerian orbital elements 

The typical planetary orbital frequencies n (see eq. (!24|) below) in the solar system 
vary from 1.3 x 10~^ Hz (Mercury) to 1.2 x 10"^'' Hz (Pluto), corresponding to timescales 
ib rangmg from 7 X 10^ s (Mercury) to 8 x 10^ s (Pluto). They are much larger than the 
time required by the light to travel across the spatial extensions of the Sun's planetary 
orbits, ranging from 2 x 10^ s (Mercury) to 2 x 10^ s (Pluto). Thus, if we consider the 
action of a monochromatic plane gravitational wave of frequency Ug and wave vector k, 
with k = Ug/c, on a planetary orbit during a time interval At comparable to an orbital 
period P^, its phase ^ = Ugt — k ■ r = Ug\t — (r/c) cos a] can be reasonably approximated by 
$ ~ Vgt, independently of the orientation a oik with respect to the planet's position r. As 
previously stated, in the rest of this Section we will assume Vg/n <^ 1 as well. 



^Notice that eq. ( 121]) has a singularity for A = 0, i.e. for a wave traveling along the 
reference x axis. 
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The standard Keplerian orbital elements of the orbit of a test particle are the 
semi-major axis a, the eccentricity e, the inclination J, the longitude of the ascending node 
Q, the argument of pericenter u, and the mean anomaly Ai. While a and e determine the 
size and the shap^, respectively, of the Keplerian ellipse, J, fi, u fix its spatial orientation. 
I is the inclination of the orbital plane to the reference {x, y} plane, while Q is an angle 
in the {x, y} plane counted from a reference x direction to the line of the nodes, which is 
the intersection of the orbital plane with the {x, y} plane. The angle u lies in the orbital 
plane: it is counted from the line of the nodes to the pericenter, which is the point of closest 
approach of the test particle to the primary. In planetary data reduction the longitude 
of the pericenter to = Q + u is customarily used: it is a "dogleg" angle. The argument 
of latitude u = w + / is an angle in the orbital plane which reckons the instantaneous 
position of the test particle along its orbit with respect to the line of the nodes: / is the 
time-dependent true anomaly. The mean anomaly is defined as 



M =n{t- tp), 



(23) 



where 



n 



^/GM/c 



(24) 



is the Keplerian mean motion related to the Keplerian orbital period by n = 27r/Pb! and 
is the time of passage at the pericenter. In the unperturbed two-body pointlike case, the 
Keplerian ellipse, characterized by 



X = r (cos Vt cos u — sin VL sin u cos /) 
y = r (sin VL cos u + cos sin u cos I) 
z = r (sin u sin J) , 



(25) 



and 



a{l-e^) 



1 + e cos / 

neither varies its shape nor its size; its orientation is fixed in space as well. 



(26) 



A small perturbing acceleration A of the Newtonian monopole, like eq. (^^, induces 
slow temporal changes of the osculating Keplerian orbital elements The Gauss equations for 



^The eccentricity e is a numerical parameter for which < e < 1 holds; e = corresponds 
to a circle. 
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their variation, valid for quite general perturbations, are (IBertotti et al.ll2003l ) 



dt 



de 
dt 

dJ 

dt 

dn 

dt 



d-uj 
dt 



nVT- 



na 



{A^sinf + Ar [cos / + i (l - ^)] } 



naVT— e 
1 



^An (^) cos 



(27) 



-Ancosf + At{1 



dM 
< dt 



n 



sin/ 



-;^^i^a)-v/i^(^ + cos/f) 



2sin^(f) 



2 ( i\ dn 

dt ' 



In eq. ([27]) p = a(l — e^) is the semi-latus rectum, and Ar, At-, A^ are the radial, transverse 
and out-of-plane components of the disturbing acceleration A, respectively. They have to 
be computed onto the unperturbed Keplerian ellipse according to 



(28) 



' An 


= A 


■R, 


< At 


= A 


•T, 


. An 


= A 


• N 



where the unit vectors along the radial, transverse and out-of-plane directions are 

cos r2 cos n — cos / sin f2 sin n, 

R= \ sin f2 cos u + cos / cos fl sin u, 

sin/ sin u, 

— cos sin n — cos / sin r2 cos n, 

T = i — sin Q sin u + cos / cos Q cos u, 

sin J cos M, 

sin/ sin 

iV = < — sin/cosr2, 

cos/. 



(29) 



(30) 



(31) 



-VI- 



In the case of eq. f|T8l) and eq. fl20l) . it turns out that it is computationally more 
convenient to use the eccentric anomaly E instead of the true anomaly /. Basically, E 
can be regarded as a parametrization of the usual polar angle Q in the orbi tal plane, bein g 
defined hA = E — esini?. To this aim, useful conversion relations are ( iBertotti et al. 



20031 ) 



cos/ 



sin/ 



dt 



cos E—e 
1— ecos_B' 



\—ecosE ' 



a(l — e cos-E), 



(32) 
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3.1. Monochromatic plane gravitational wave propagating along a generic 

direction 

Straightforward calculationl^ yield the long-term, i.e. averaged over one orbital period, 
variations of all the Keplerian orbital elements of the test particle due to eq. ([7]). They are 

^ = 0, 

5e>/l-e2 

- TFti { ~^ ^ '^^^ {^xs COS f2 + /Cyz sin H) + 

+ 4 cos / cos 2a; [ — 2^xy cos 2n+ {fC^x ~ K,yy) sin 2flj + 

+ sin 2u! ^(/Cxa; — l^yy) {3 + cos 2/) cos 2f2 + 6 (/Cxx + l^yy) sin^ /+ 

+ 4 sin 2/ (Kjci sin Q - Ky^ cos fi) + 2)Cxy (3 + cos 21) sin 2n] } , 

^ = -^^ - [|4cos / (^2 + 3e2 + Se^ cos 2a;) (Kxz cos Q + K:,,^ Sinn) - 

- 2 (2 + Se^) sin / [2Kxy cos 2n + (Ky^ - K^x) sin 2n] - 

- lOe^sin/ cos 2a; [2K^xyCos2n + {^yy ~ K^xx)sin2n] + 
+ he^ sin 2a; |4 cos 2/ {K,yz cos n — ^xs sin H) + 

+ sin 2/ [-3 (Kxx + ICyy) + (Kxx - ICyy) COS 2n + 2Kxy sin 2n] } I] , 

^ = -^1 = |2cos2/csc/ (^-2 - 3e2 + Se^ cos2aj) (Kyj cosn - Kxj Sinn) + 

+ cos / (-2 - 3e^ + 5e^ cos 2a;) [-3 (Kxx + Kyy) + (Kxx - Kyy) cos 2n + 2K,xy sin 2n] + 

(33) 

+ he^ sin 2a; [2K^xy cos 2n — 2 cot / (/Cxs cos H + K^ys sin H) + {ICyy — K^xx ) sin 2n] j- , 
^ = -L = |-3 (Kxx + Kyy) (1 - 4cos / + cos 2/) (^-2 - 3e^ + 5e^ cos 2a;) + 

+ (cos / - 2) cos2n [2 (Kxx - Kyy) cos / (-2 - 3e2 + Se^ cos 2a;) + 20e2K:xy sin 2a;j + 

+ 2 (cos / - 2) cos n CSC / [2K:yi cos 2/ ^-2 - 36^ + Se^ cos 2a;) - lOe^ Kxj cos / sin 2aij + 

+ 2 (cos 1-2) Sinn esc / [-2Kxz cos 21 (-2 - 3e^ + 5e^ cos 2a;) - lOe^Kyj cos / sin 2a;j + 

+ 2 (cos / - 2) sin2n [2Kxy cos / (-2 - 3e^ + 5e^ cos 2a;) + 5e^ (Kyy - Kxx) sin 2a;j } , 
= „+ |8 [-7 - 3e2 + 5 (1 + e^) cos2a;j sin2/ (Kyj cosn - Kxz Sinn) - 

- 80 (1 + e^) sin /sin 2a; (Kxj cos n + Ky^ sin n) + 

+ 2 cos 2/ [-7 - Se^ + 5 (l + e^) cos 2a;j [3 (Kxx + ICyy) + {ICyy - Kxx) cos 2n - 2K:xy sin 2n] - 

- 2 1^7 + 3e^ + 5 (1 + e^) cos 2a;j [Kxx + Kyy + (Kxx - Kyy) cos 2n + 2K:xy sin 2n] - 

- 40 (1 + e^) cos / sin 2a; [2K:xy cos 2n + (Kyy - Kxx) sin 2n] } . 

In order to deal with manageable expressions of general validity, we did not display in 
eq. f l33p the explicit expressions of the tidal matrix coefficients in terms of /3 and A; The 
result of eq. f l33p is exact both in e and in I in the sense that no a priori simplifying 
assumptions about the eccentricity and the inclination of the perturbed test particle's orbit 
were assumed in the calculation. It can be noticed that the semi-major axis is not affected 



^We do not explicitly show here the analytical expressions for Ar, At, An due to eq. f|T8|) 
because they are rather cumbersome. 
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by the passage of a very slowly varying gravitational wave; the variation of the eccentricity 
is proportional to e itself, so that a circular orbit does not change its shape. Moreover, 
eq. ( I33l) does not contain secular effects because of the presence of /, Q, u which, actually, 
experience slow change^ in time as far as the planets of the solar system are concerned. A 
further modulation is due to the coefficients ]Cij of the tidal matrix containing the (slow) 
harmonic variation of the impinging gravitational wave. Since all such frequencies are much 
smaller than the typical orbital ones for the Sun's planet, the terms containing them were 
kept fixed in the integration yielding eq. 



3.2. Monochromatic plane gravitational wave propagating along a the z axis 

By using the same procedure it is possible to work out the long-term variations of 
the Keplerian orbital elements for a direction of incidence of the wave coinciding with the 
reference z axis. By using eq. (120|) one gets 



r da 

dt 



0, 



de cos / COS 2uj{h2 cos 2Q— hi sin 2f2)+2(3+cos 21) sin 2aj(hi cos 2fl+h2 sin2f2)] 

H ~ 16n ' 

sin/|^— (2+3e'^+5e^ cos 2a;) {h2 cos 2fl—hi sin2n)+5e^ cos /sin 2a; ( ft, i cos 2n+ft2 sin2f2)] 
'dt ~ inVl-e'^ ' 

5e^ sin2a;(ft2 cos 20— fti sin 2f2)+cos /(— 2— 3e-^+5e^ cos 2a;) (hi cos 2fl+h2 sin 2f2) 
'dt ~ 4nVl-e^ ' 

^ = Ji„,. {-40 [e^ + (e^ - 2) cos /] sin 2uj (/la cos 2Q - hi sin 2Q) + , ^^^^ 
+ 20 (3 - 2e^ - 2e^ cos / + cos 21) cos 2uj {hi cos 2fi + /12 sin 2fi) + 
+ 16(3-3e2 + 5cos/)sin2(|) (/ii cos 2Q + /i2 sin 2Q) } , 
= n+ 3^[-80(l + e2)cos/sin2a;(/^2Cos2fi-/^lsin2fi)- 
- 20 (1 + e^) (3 + cos 21) cos 2lo {hi cos 2fi + /12 sin 2fi) - 
- 8(7 + 3e2) sin^ I {hi cos 2fi + /i2 sin 2fi)l . 



Also in this case, no a priori assumptions on e and / were made, so that the rates of eq. 
( 134|) are exact in this respect. Notice that for / = 0, i.e. for incidence of the gravitational 



^''They are mainly due to the mutual N— body classical perturbations: the rates of change 
can be found at jhttp: //ssd.jpl. nasa.gov/txt7p_elem_tl.txt 
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wave normal to the orbital plane, the inclination is left unaffected. Moreover, also in 
this case, the elements hi and h2 of the tidal matrix, which are (slowly) time-varying 
harmonic functions whose amplitudes are proportional to z/|, were assumed constant in the 
integrations over the planet's orbital period. 

Notice that, for a plane wave traveling along the z axis, i.e. for eq. (fT9|) . eq. (!33|) 
reduces just to eq. It is straightforward to obtain the formulas valid also for the other 

specific directions of propagation examined in Section |5] by suitably specializing eq. (15^ to 
such cases (cfr. eq. ( IT2l) . eq. ( IT3l) and eq. (fT5l) ). 



3.3. A comment on the approximation used for the harmonic wave functions 

In obtaining eq. (!33|) and eq. ( 134|) we kept the tidal matrix coefficients, which include 
the time-dependent harmonic functions of the wave, constant in the integrations over 
one orbital revolution of the test particle. This implies that we considered waves having 
frequencies Ug much smaller than the orbital ones n. As we will see in Section HI it is 
in contrast with the other approaches followed so far in literature, in which no similar 
assumptions on Ug were made. Apart from the fact that, from a computational point of 
view, our approach allows for exact calculations in e and J, our choice is justified also 
from a phenomenological point of view. Indeed, applying our results to the solar system 
implies that we could, in principle, put constraints over a part of the ultra-low frequency 
spectrum of gravitational waves for which neither ground-based nor space-based dedicated 
experimental devices are availablj"]. Moreover, in view of continuous tracking of planets 
by means of ranging either directly to their surfaces or to ever more numerous orbiting 
spacecraft it will be possible, in principle, to obtain more and more accurate bounds 
because of increasing data records over the years. On the other hand, it is, after all, of little 
practical utility to perform cumbersome calculations of the wave-induced orbital effects 
involving frequencies as large as, or even larger than, the planet's ones since much more 
accurate dedicated experiments already exist covering such windows of the spectrum of 
gravitational wav es. Indeed, the accuracy reachable in future p l anned interplanetary la ser 



ranging facilities (jSmith et al.ll2006t iTuryshev fc Williamsl 120071 : iMerkowitz et al.l 120071 ). of 



the order of about 1 — 10 cm, is not comparable with that of the latest gravitational wave 
detectors; it can be acceptable in order to put upper bounds when no other, more accurate 
devices exist, as for the ultra- low frequency waves. 



-^^A possible exception may be the proposed Inflation Probe by NASA (iBenford et al. 
20091 ). dedicate d to map the polari zation of the Cosmic Microwave Background at 10"^^ Hz 



(see Figure 1 of iPrince et al.l ( l2010l )). Notice also that precision timing of millisecond pulsars 
may be used for gravit ational waves in the range 10~^ — 10~^ Hz. (IKopeikinl 119971 : iHand 



20101 : 1.Tenet et al.ll201lh . 
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4. Confrontation with other approaches in hterature 



In this Section we briefly outline the approaches followed by some other researchers, 
with particular care to the orbital effects caused by the traveling wave. 



Rudenkd (119751 ) uses polar coordinates in the orbital plane and considers an orbit 
disturbed by a monochromatic plane gravitational wave traveling along a generic direction; 
the additional force resulting on the test particle is obtained within the linear approximation 
of general relativity (within the frame work of t h e "th eory of gravitation in plane space" by 
Zel'dovich fc Novikovl ( 1l967l )). Then, iRudenkd (119751 ) assumes a normally incident wave 
and solves the equations of motion at zero order in eccentricity. Finally, he studies the 
variation of the orbital radius Ar for various values of the wave's frequency. 



Turnerl ( 119791 ) works in the TT gauge by deriving the geodesic equations of motion 
in cartesian coordinates of both the binary system's constituents. Then, he takes their 
difference obtaining the equations for their relative motion. It turns out that th ey are 
formally different from those obt ained by other authors like, e.g.. iMashhooru (119781 ) from the 



geodesic deviation equation, but iTurnerl (jl979[ ) shows th at they are, ac tually, equal up to an 
unphysical coordinate transformation. In the scenario by iTurnerl (11979! ) . the monochromatic 
plane wave travels along the negative z direction. He uses the Gaussian perturbative scheme 
to work out the temporal changes over one orbital revolution of the semi-major axis and the 
eccentricity of a circular orbit by settin g ^ = w = / = 0. Then, he discards the limitation 
/ = 0, but not the other ones. Finally, iTurnerl ( 1l979l ) considers non-circular orbits, and 
computes the variation of a to the lowest order in e for f2 = / = = 0. In all of such cases, 
he takes different values for Vgjn. 



Mashhoon et al.l (jl98ll ). relying upon the methods developed in iMashhoonI (jl978[ ). look 
at the effects that various kinds of gravitational radiation, among which monochromatic 
plane waves are c onsidered as well, have on a Newtonian binary system. Let us recall that 
MashhoonI (119781 ). in considering slow particle motions and weak waves having wavelengths 
larger than the syst em's orbita l size, adopts the linearized Jacobi equation for the geodesic 
deviation equation. iMashhoonI (119781 ) writes down the equations of motion in cylindrical 
coordinates and solve them for various values of Vg with respect to n. A s far as the 
monochromatic plane wave case is concerned, also iMashhoon et al.l (jl98ll ) assume that its 
wavelength is much larger than the size o f the binary syste r n; no simplifying assumptions 
on the wave's frequency are made. Then, iMashhoon et al.l (Il98ll ) work out the resulting 
relative change Ar/r occurring in the test particle's distan c e from the primary in the 
low-eccentricity approximation. Moreover, iMashhoon et al.i (Il98li ) consider also the case 
of a pair of planets moving in the same plane along circular orbits, and calculate the 
wave-induced relative change ATZ/TZ in their mutual distance. 



Nelson &: Chaui (119821 ) look in the TT gauge at a monochromatic plane gravitational 



wave traveling along the z axis, and having wavelength much larger than the size of the 
system considered. They write down the components of the wave-induced acceleration 



-17- 



experienced by the test particle in cylindrical coordinates for generic shape and inclination of 
the orbit. Then, they specialize the m to the case of normall y incident wave, i.e. for an orbit 
with zero inclination. At this point, iNelson &: Chad (119821 ) solve the resulting perturbed 
equations of motion by evaluating the radial change Ar over one orbital revolution. In their 
computation, they resort to some approximations in e, and do not consider the harmonic 
wave functions constant over typical orbital timescales. Th e resulting changes for different 
values of Vg/n are inspected. Then, iNelson fc Chaul (119821 ) examine the case of a generic 
inclination of the test particle' s orbit by worki n g out the shift along the wave's propagation, 
i.e. along the z axis. Finally, INelson &: Chaul ( 119821 ) perform some numerical analyses of 
the radial shifts for different values of the eccentricity and of the phases of the wave. 



Also llvashchenkol ( 119871 ) considers the case of a monochromatic plane gravitational 
wave, in the TT gauge, traveling along the z axis. He assumes that its wavelength is much 
larger than the characteristic size of the pertur bed system, so that he can neglect the term 
A; • r in the wave's phase. On the other hand, llvashchenkol (119871 ) does not make any a 
priori assumptions about the magnitude of the wave's frequency Ug with respect to the 
orbital one n. After having obtained the exact radial, transverse and no rmal components o f 
the perturbing acceleration from its expression in cartesian coordinates, llvashchenkol (119871 ) 
makes use of the Gauss equations for the variation of all the Keplerian orbital elements, 
apart from the mean anomaly. At this point, after having set w = 0, he makes an expansion 
of the right-hand sides of the Gauss equations, evaluated onto the unperturbed Keplerian 
ellipse, to first o rder in e by using t he mean anomaly M. as independent variable. Thus, in 



view of eq. (|23|) , llvashchenkol (119871 ) obtains approximate expressions of the form 



dt 



d{t), ^ = a,e,I,n,u, 



(35) 



which he straightforwardly integrates. In doing that, he does not con sider and h2 as 
constant by keeping, instead, their harmonic time dependence. Then, llvashchenkol (119871 ) 
discusses various particular cases for different orbital geometries and, especially, for different 
values of i^g/n. 



Also iKochkin fc Sbytovl ( 119871 ) adopt a wave in TT gauge traveling along the reference z 
axis. Moreover, they assume for it a circular pol arization, and do no t make any assumptions 
about the wave's frequency and its wavelength. iKochkin fc Sbytovl (Il987l ) use cylindrical 
coordinates to write down the spacetime line element for a linearized superposition of the 
Schwarzschild and wave fields, from which the particle's equations of motion are obtained. 
Concerning the orbital geometry, they put the perturbed orbit , not assumed a pr i ori 
circular, in the reference {x, y} plane, so that k • r = 0. Then, IKochkin fc Sbytovl (119871 ) 
do not rec ur to any kno\vn per turbative schemes, like, e.g., the Gauss and the Lagrange 
equations ( iBertotti et al.ll2003l ). in dealing with the resulting equations of motion. After 
some changes of variables, they integrate the resulting modified equations over one orbital 
revolution by inferring the cumulative changes of the semi-major axis, th e eccentricity and 
the pericenter after one orbital period. Finally, IKochkin fc Sbytovl (Il987l ) consider how the 
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results they obtained vary for different values of z/^/n; they consider the circular case as 
well. 



Chicone et al.l (Il996al ) and IChi cone et al.l f ll996bl ) look at the action of a monochromatic 
plane gravitational wave traveling along the reference z axis and normally impinging on 
an unperturbed two-body system. The wave's wavelength is assumed to be much larger 
than the semi-major axis of the orbit, and the characteristic velocities within the system 
are much smaller than c. No assumptions are made a priori on the wave's frequency Vg. 
The resulting relative accele ration is obtained from the geodesic deviation eq uation in 
cartesian coordinates. Then, IChicone et al.l (jl996al ) and IChicone et al.l ( 1l996bl ) write down 
the Hamiltonian of the perturbed system in polar coordinates, and adopt the Delaunay 
orbital elements L, G, £, g for a non-circular orbit. Fourier series expansions in terms of the 
mean anomaly are performed. Different values for Vgjn and wave's polarizations are, then, 
considered. 



Chicone fc MashhoonI ( 12002| ). after having generalized the Jacobi equation by taking 
into account also the relative velocity of the geodesies, consider a monochromatic plane 
gravita tional wave traveling along t he reference x axis of the local Fermi quasi-inertial 
frame. IChicone fc MashhoonI ( 12002| ) write down the resulting equations of motion of the 
test particle in cartesian coordinates; they are r ather involved because of the presence of 
the non-linear, velocity-dependent terms. Then, IChicone fc MashhoonI ( 120021 ) study the 
motion of the test particle along the x axis itself. 



The general wave-binary s ystern scenario adopted by llsmaiel fc SaadI (120111 ) is analogous 



to that by, e.g., Ilvashchenkol (119871 ). However, after havi ng written the compon ents of 



the tidal-type wave acceleration in cartesian coordinates, llsmaiel fc SaadI (120111 ) obtain 
a potential function 7?. from t hem in order to use the perturbative scheme based on the 
Lagrange planetary equations (IBertotti et al.ll2003l ) for all the Keplerian orbital elements, 
apart from M.. Then, they evaluate 7^ onto the unperturbed Keplerian ellipse, without 
fixing a;, and expand it to some, unspecified, order in e by using the mean anomaly a 
independent variable. In such a way, they obtain that the right-h and sides of the Lagrange 
equations depend only on t, which allows llsmaiel fc SaadI (120 111 ) to integrate them over 
one orbital period. Notice that they do not keep the wave's harmonic function constant in 



the integration. Finally, they compute the long-term changes of I a, e, /, f2, uj of Venus and 
Pluto for three, unspecified, different sources of gravitational waves. 



^^It is unclear if the figures in Table III of llsmaiel fc SaadI ( 120111 ) refer to the rates of 
change of the Keplerian orbital elements, i.e. if they are m s~^ and deg s~^, or to their shifts 
over one orbita l perio d, i.e. if they are m and deg. Also the units adopted in Table III of 



Ismaiel &: SaadI (120111 ) are unspecified. 
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5. Summary and conclusions 

We analytically worked out the long-term variations of all the six Keplerian orbital 
elements of a test particle orbiting a central body due to the perturbing tidal acceleration 
induced by the passage through the system of a monochromatic plane gravitational wave. 
We assumed that its frequency is much smaller than the orbital frequency of the test 
particle, so that the time-dependent harmonic functions of the wave were kept constant in 
the integrations performed over one orbital period. We considered a generic direction of the 
wavevector. 

We obtained that all the Keplerian orbital elements, apart from the semi-major axis, 
undergo long-term variations. They are linear combinations of the non-vanishing elements 
of the tidal matrix, containing the slowly time-varying harmonic terms due to the wave, 
with coefficients which are complicated functions of the eccentricity, the inclination, the 
node and the pericenter of the test particle. We did not make any a priori simplifying 
assumptions concerning the inclination and the eccentricity of the test particle's orbit, thus 
obtaining exact results as far as such a point is concerned. The variation of the eccentricity 
is proportional to the eccentricity itself, so that circular orbits do not change their shape. 
In the case of incidence normal to the orbital plane, also the inclination is left unaffected. 

From a practical point of view, in the most general case one has six unknowns: the 
two angles of the wavevector, the frequency of the wave, its two independent amplitudes, 
and a phase lag. In principle, it is possible to constrain all of them by determining 
the secular variations of some Keplerian orbital elements for, e.g., several planets of the 
solar system. Actually, this seems to be the current trend in astronomical research since 
extensive planetary data reductions, in which the corrections to the standard precessions of 
the perihelia and the nodes of some major bodies of the solar system are estimated, are 
currently ongoing by independent teams of astronomers. Future implementation of planned 
or proposed accurate interplanetary ranging facilities may further enhance such a program. 
It should also be considered that no other dedicated experiments are currently available for 
the ultra-low frequency part of the spectrum of gravitational waves. 

Our results are quite general since they are valid not only for the solar system's planets, 
but also for any generic gravitationally bound two-body systems whose characteristic 
orbital frequencies are quite larger than the incident wave's frequency. Moreover, they can 
be extended to any generic perturbing acceleration of tidal type with constant coefficients 
over the characteristic timescales of the perturbed system. 
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